Renormalized spectral function for Co adatom on the Pt(lll) surface 
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The strong Coulomb correlations effects in the electronic structure of magnetic Co adatom on 
the Pt(lll) surface have been investigated. Using a realistic five d-orbital impurity Anderson 
model at low temperatures with parameters determined from first-principles calculations we found 
a striking change of the electronic structure in comparison with the LDA results. The spectral 
function calculated with full rotationally invariant Coulomb interaction is in good agreement with the 
quasiparticle region of the STM conductance spectrum. Using the calculated spin-spin correlation 
functions we have analyzed the formation of the magnetic moments of the Co impurity orbit als. 

PACS numbers: 71.27.+a,73.20.At, 75.20.Hr 
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Nanosystems consisting of magnetic atoms deposited 
on nonmagnetic surface become attractive for the devel- 
opment of novel types of memory and computation de- 
vices. The electronic and magnetic properties of such 
systems can be probed by means of the inelastic scan- 
ning tunneling spectroscopy (STM) [1]. The STM spec- 
tra demonstrate different types of elementary excitations. 
For instance, if the magnetic atoms are deposited on in- 
sulating surface there are steps structure in the differen- 
tial conductance dl/dV spectrum [2] that is related with 
inelastic excitations and can be reproduced in the frame- 
work of the spin tunneling Hamiltonian approach [J, |4[ . 
The latter is based on the fact that well-localized mag- 
netic moments of the adatoms can be described with a 
simple Heisenberg spin-Hamiltonian 0, [H[ . 

In contrast to the insulating substrate, the experimen- 
tal STM spectrum of magnetic adatoms deposited on a 
metallic surface can show the Kondo-like physics @, . 
Nevertheless, the recent spin-polarized STM experiments 
[8, 9] for Co impurities on Pt(lll) surface have demon- 
strated the existence of narrow peaks near the Fermi level 
which are not split in the external magnetic field. More- 
over, the width of this peak of around 150 meV is much 
larger than that one would expect for a Kondo system 
[si EqI • From the theoretical point of view this peak can 
be related with the local density of state of the magnetic 
adatom 11]. Indeed, the DFT-based calculations lJ-Q 
have demonstrated a peak near the Fermi level, however, 
the width of the peak is much larger than the experi- 
mental one and is located fare away from zero bias. The 
main problem with the DFT-calculations is the missing 
of dynamical quantum fluctuations which can be very im- 
portant in case of strongly correlated nano-systems [l5| . 
They can lead to a considerable renormalization of the 
LDA spectrum near the Fermi level [16| . In order to cap- 
ture the local correlations on magnetic adatom properly 
one should solve the many-body Anderson-type Hamil- 
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FIG. 1: (color online). Densities of states obtained from LDA 
calculation (brown line) and many-body Anderson model cal- 
culations (bold blue line) in comparison with of the quasipar- 
ticle region of the experimental differential conductance 0] 
(thin gray line). 



tonians @, Ell at experimental temperatures. 

In this Letter we solve the realistic five d-orbital An- 
derson impurity model with parameters determined from 
first-principles calculations to describe the electronic and 
magnetic properties of Co/Pt(lll) system. Using differ- 
ent parametrization of the Coulomb interaction matrix 
we have analyzed the microscopic origin of low-energy 
excitations of the system. Our main result can be in 
short presented in Fig.l which gives the comparison be- 
tween the quasiparticle region of the experimental STM 
differential conductance [8[ with the spectral functions 
obtained from the LDA and temperature-dependent An- 
derson impurity-model calculations. One can immedi- 
ately realized that there is a huge renormalization of 
the quasiparticle density of states in Anderson model 
in comparison with the LDA results. Moreover a good 
agreement between experimental results and Anderson 
model can be obtained only with the full rotationally- 
invariant parametrization of Coulomb interaction matrix. 
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The peaks in the density of states below and above the 
Fermi level are formed by Co-3d states of different sym- 
metry which can be checked in the experiments. We have 
also calculated the partial spin-spin correlation functions 
for different Co orbitals in order to study formation of 
the magnetic moment. 

LDA analysis. We performed the LDA calculation 
with structural relaxation for single Co atom on the 
Pt(lll) surface by means of the projected augmented- 
wave (PAW) method [18] as implemented in the Vienna 
ab-initio simulation package (VASP) [19j. We used an 
energy cutoff of 300 eV in the plane-wave basis con- 
struction and the energy convergence criteria of 10 ~ 7 eV. 
Atomic positions of the considered system were relaxed 
with residual forces less than 0.01 eV/A To simulate the 
atomic structure of Co adatom on the Pt(lll) surface we 
have used the supercell approach. The supercell contains 
three-layer (3 x 3) Pt(lll) surface, Co atom and vacuum 
region of 10 A. Lattice constant for Pt lattice was cho- 
sen to be 3.92 A that is the experimental value of lattice 
constant for the bulk fee Pt. The lowermost layer of Pt 
has been fixed under relaxation. The obtained vertical 
distance of Co atom of 1.82 A is in agreement with the 
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Fig. 2 gives the LDA partial densities of states. All the 
cobalt 3d states can be classified with respect to the hy- 
bridization with Pt states. One can see that the density 
of states of 3z 2 — r 2 symmetry lies below the Fermi level 
and demonstrates a localized peak. It is due to a small 
overlap and hybridization with states of the nearest Pt 
atoms. In turn the states of the in-plane orbitals of xy 
and x 2 — y 2 symmetry hybridize with Pt states and are 
wider than that of 3z 2 — r 2 orbital. The largest hybridiza- 
tion occurs between xz (yz) and 5d Pt states. 

To study the magnetic properties of Co/Pt(lll) we 
have performed spin-polarized LDA calculations. The 
obtained value of the local spin magnetic moment Ms = 
2.09 jiB is the same for out-of-plane and in-plane direc- 
tion of magnetization, whereas the orbital moment varies 
from M[ = 0.16 /i^ to M£ = 0.09 These values are 
much smaller than the total giant magnetic moment of 
Co-Pt systems observed experimentally. 

Coulomb correlations. One can see that in the frame- 
work of the LDA the orbital polarization is strongly un- 
derestimated. The same result was obtained in the pre- 
vious first-principles investigations [l3|, QjJ . To overcome 
this problem the different kinds of the orbital polarization 
corrections were used to fit the experimental data. Us- 
ing the static Coulomb correlations within the LSDA+U 
method the correct value of the orbital magnetic moment 
in Co-Pt systems have been reproduced [12 1. 

The static mean-field LDA+U and many-body ap- 
proaches we use require the on-site Coulomb interaction 
and intra-atomic exchange interaction as input param- 
eters. To calculate them we have used the constrained 
LDA procedure [20|. Our calculations have shown that 
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FIG. 2: (color online). Densities of states obtained from LDA 
calculations (gray solid line). The blue dashed lines corre- 
spond to fitting results obtained by using Eq.(2). 



the values do not differ significantly for relaxed (U = 
6.6 eV, J = 0.9 eV) and non-relaxed (U = 6.9 eV, J = 
0.9 eV) structures. Thus in our many-body calculations 
we have used the averaged values U= 6.75 eV and J = 
0.9 eV. Using the estimated U and J we have performed 
the LDA+U calculations in the framework of the tight- 
binding linear muffin-tin orbital approach within the 
atomic sphere approximation (TB-LMTO-ASA) method 
2lL |22| . The calculated value of the magnetic moment 
Ms = 2.8 /ib is larger than those experimentally ob- 
served [23(1 . It is well known that the account of the 
many-body correlations can result in a renormalization 
of the electronic spectrum and as a consequence the mag- 
netic moments becomes smaller [16[ . 

Many-body calculations. The importance of account of 
the dynamical Coulomb correlations for description of the 
adatom on a metallic surface was emphasized in Refll7. 
The authors have investigated the electronic structure 
of a cobalt impurity on the surface of Cu(lll) by using 
a CT-QMC method. They have found that in contrast 
to the bulk system there is a large sign problem related 
with non-diagonal elements of the Coulomb interaction 
matrix. The sign problem prevented them from calculat- 
ing the physical properties at low temperatures. 

Here we use the finite temperature exact diagonal- 
ization Lanczos-method [24( in order to solve the An- 
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derson impurity Hamiltonian for realistic multi-orbital 

systems 
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25j . We consider such a method to be prefer- 
able in comparison with others because one can use the 
full U-matrix and the simulation can be performed at ex- 
perimental temperatures [23| . The main disadvantage of 
the exact diagonalization is a discretization of the LDA 
spectrum required to define the impurity and bath ener- 
gies as well as the hopping parameters in the Anderson 
Hamiltonian. 

To study the correlation effects we diagonalized the 
Anderson impurity Hamiltonian [26| 
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where (e/~) is the energy of the impurity (surface) 
states, <i+ a (c^ a ) is the creation operator for impurity 
(surface) electrons, Vmk is a hopping between impurity 
and surface states, and U ram 'm"m"' is the Coulomb ma- 
trix element [17j . Importantly, in our calculations the full 
Coulomb interaction vertex Eq.(l) as well as the density- 
density form H den \Y,mm'aVmm'mm'n Tna n Tn >- (7 + 

(Umm>mm> U mm > m > m )n ma n m > a } were used. The chemi- 
cal potential \i was fitted to obtain the number of impu- 
rity electrons of 7.6. 

The energies e m and the hopping parameters V m k 
were found by minimizing the LDA Green functions using 
the following expression 
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In our calculations we used N s — 2 for each impurity or- 
bital, thus the total number of the electronic levels was 
equal to 15. The same number of the bath and impurity 
orbitals was used in Ref. [2?] to solve a DMFT prob- 
lem. Fig. 2 gives the comparison between LDA densities 
of states and fitted spectral functions. One can see that 
there is good agreement between them. The obtained 
values of e m agree well with an energy centers of the cor- 
responding orbital. The smallest hybridization hoppings 
were obtained for 3z 2 — r 2 orbital, since this orbital is the 
most localized one. 

The impurity Green's functions and spin susceptibili- 
ties are of our interest and we calculate them by means 
of the general expression for the correlation function in 
the Lehmann's spectral representation: 
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where X™ n , is a matrix element of the fermion (£ = 1) 
or spin (^=-1) operator and Z is the partition function. 
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FIG. 3: (color online). Densities of states obtained from 
model calculations for U = 6.75 eV and J = 0.9 eV. The 
dashed gray and solid blue lines correspond to calculations 
with density-density Coulomb interaction and full U-matrix. 
The Fermi level is zero. 



To avoid the extrapolation from Matsubara frequencies 
to real frequencies we use energy lj + i5, where 5 = 
At low temperatures the small number of excited states 
contribute to the Boltzmann factor in Eq.(3). In our in- 
vestigation we calculated these correlation functions for 
40 excited states, which is enough to work in the temper- 
ature range - 50 K with the Boltzmann factor smaller 
than 10~ 6 . 

In order to calculate the impurity Green's function 
(single-particle correlation function) we used X™ n , = d nn > 
in Eq.(3). These results for different types of Coulomb 
interaction matrices are presented in Fig. 3. The cal- 
culations with density-density and full Coulomb interac- 
tion matrices give slightly different values of the magnetic 
moment of 2.2 \±b and 2.0 respectively. It means 
that the account of the non-diagonal elements of the U- 
matrix results in an additional suppression of the mag- 
netic state. One can see that in the case of the density- 
density Coulomb matrix the spectrum near the Fermi 
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FIG. 4: (color online). Impurity z-z (left panel) and total 
spin-spin (right panel) correlation functions obtained from 
Anderson model calculations with full Coulomb interaction 
matrix. Solid gray, dashed light blue and dotted red lines 
correspond to spin susceptibilities of 3z 2 — r 2 , xy (x 2 — y 2 ) 
and yz (xz) orbitals, respectively. 

level is strongly renormalized due the spectral density 
transfer to the empty states at + U. In turn the ac- 
count of non-diagonal elements of the U matrix results 
in new excitations of the system near the Fermi level. 

The density of states obtained from the Anderson 
model calculations with the full U-matrix is in good 
agreement with the quasiparticle region of the differential 
conductance dl/dV (Fig.l). There is a peak at about - 
0.05 eV below the Fermi level. This peak results from xy, 
yz, xz and x 2 — y 2 states which are strongly hybridized 
with conductance states of the Pt(lll) substrate. The 
width of the peak is 100 meV, which is much larger than 
that one obtained for a Kondo system. The states just 
above the Fermi level are mainly of 3z 2 — r 2 symmetry 
that is hybridized much less than other. Such a sym- 
metry difference of the states near the Fermi level can 
explain the different behaviour of the conductance at the 
external magnetic field Q. When the magnetic field is 
reversed the peak below the Fermi level strongly changes 
the intensity. It is due to the fact that the peak is origi- 
nated from the partially occupied cobalt states that give 
the largest contribution to the total magnetic moment. 

Another important source of information about an im- 
purity system is the spin susceptibility that can be cal- 
culated by using Eq.(3) with X™ n , = S* n , and £=-1. 
These results are presented in Fig. 4 (left panel). One 
can see that all the susceptibilities yield peak which is 
mostly pronounced for xy(x 2 — y 2 ) and 3z 2 — r 2 states. It 
means that there is the local moment for these states [28|. 
As for the yz(xz) state the susceptibility is nearly flat 
and we believe that the spin correlation function demon- 
strates an itinerant behaviour due to a strong hybridiza- 
tion between Co and Pt states. We have also calcu- 
lated the correlation function for the linear combination 

X™ n 'X™ n = S* n ,S*r n + \(Snn'Sn'n + ^nn'^n'n)' ^he 

resulting correlation functions presented in Fig. 4 (right 
panel) demonstrate high-energy excitations which are 
mainly of 3z 2 — r 2 symmetry. One can use the calcu- 
lated spin susceptibilities in order to simulate inelastic 



excitations which are related with the transition matrix 
elements for the spin operators [H, 0, Eq| . We left such 
an analysis for a future study. 

In conclusion, the investigation of the electronic and 
magnetic properties of Co/Pt(lll) system was carried 
out in the framework of the impurity Anderson model 
by using finite temperature Lanczos-solver. Using the 
obtained results we proposed an interpretation of the ex- 
perimentally observed energy spectrum of the magnetic 
nanosystem on the surface in terms of strong many-body 
correlation effects. 
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